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Abstract 



We present study of finite-size scaling and universality of crossing probabil- 
ities for the (/-state Potts model. Crossing probabilities of the Potts model 
are similar ones in percolation problem. We numerically investigated scaling 
of ir s - the probability of a system to percolate only in one direction for two- 
dimensional site percolation, the Ising model, and the q-state Potts model for 
q = 3, 4, 5, 6, 8, 10. We found the thermal scaling index y = - for q < 4. In 
contrast, y ^ - for q = 4. 

I. INTRODUCTION. 

In recent years the scaling and universality of the g-state Potts model and percolation are 
subject of intensive study pHJ]. There are two kinds of quantities of interest in percolation 
problem: cluster distribution functions and crossing (spanning) probabilities. With n s (p) - 
the number of clusters of size s per lattice site, the m{p) = s max n Smax (p) = p(l— sn s (p)) 

is probability of a lattice site to belong to "infinite" cluster (the cluster of maximum size) 
and x{p) = X] s 2 n s (p) is the mean size of a finite cluster. For an infinite lattice near 

S 7^ S m a x 

the percolation point p c xip) ~ \p ~ Pc|~ 7 , m ~ (p — Pc) 13 , P > Pc @- This quantities 
corresponds to magnetization m{T) and magnetic susceptibility x(T) of the Potts model. 

Critical properties of percolation may be found from the study of crossing probabilities 
7r(p; L). Crossing probabilities were used by Reynolds, Stanley and Klein for renormalization 
group study of percolation J/]J§. Universality of this crossing probabilities as functions of 
an aspect ratio r in the critical point was found in |J. 

The exact results for ?r^(r) and n hv (r) in the percolation point was developed by 



Cardy [J^] and Watts |TT[ for different r. 

Scaling properties of the function nh(p',L) - the probability, that a system percolates 
in the horizontal direction for two- and three-dimensional lattices was investigated by Hu, 



Lin, Chen |I2-LH]. They show, that ~Kh{p) is a universal function of scaling argument x = 
(p — Pc)L» , where p c - critical point, L - lattice size, v - correlation length exponent. 

Universality of the crossing probability 7r^,(/3) for Ising model was investigated by Lang- 
lands et. all. [|l^ 



In this paper we present results of numerical investigations of the crossing probability 
7r s (/3; L) (the probability of a system to percolate only in the one, horizontal either vertical, 



1 



direction) for the g-state Potts model. We found, that ir s (/3;L) is a universal function of 
the scaling variable {(3 — /3 c )L y , where for the Ising model and the q = 3 Potts model the 
thermal scaling power y = -. Further, the limit of the crossing probability at the critical 
point tt s {Pc'i L) is nonvanishing for the Ising model and for the Potts model with q = 3, 4. In 
the case of the Potts model q — 5,6, 8, 10 percolation probabilities goes to zero, when the 
system size L tends to infinity. 

The crossing probability in the critical point behaves as tt s (/3 c ;L) ~ A — alog(L) for 
L < £ c , while the dependence is exponential 7T s (f3 c ; L) ~ Aexp(-aL) for L > £ c , where £ c is 
the correlation length in the critical point. 



II. CROSSING PROBABILITIES. 

A. Crossing probabilities for percolation. 

Let us consider the site percolation on the square lattice. Each site is occupied with 
probability p and is empty with probability 1 — p. Let us denote by u the sample, i.e. the 
fixed distribution of occupied sites on the lattice. The full set of samples fl consists of 2 L 
configurations. 

Let us introduce indicator functions of crossing Ih{w) in horizontal and I v {uS) in vertical 
directions For example, if the cluster spans the sample horizontally then Ih(u) = 1 

otherwise Ihioj) = 0. Let us S(u) - the number of occupied sites in the configuration u>. 

The crossing probabilities for the site percolation on the square lattice are introduced by 

n h (p;L) ^P*H4H = Xy M (l -p) L2 - S ^I h (cu) (1) 
7r v (p;L) = Y^PsiteWM = Xy (h °(l -p) £a - 5H /„(w) (2) 
7r hv (p;L) = Y,PsitMhW v {u) = YjP^^-vf'-^hlM (3) 

vr s (p; L) = J2p SH (1-p) L2 - S(uj) - /.(w)) + (1 - hW v {oo)\ (4) 

wen 

The function 7Th(p] L), defined by expression ([l|), is the probability, that the lattice with 
linear size L percolates in the horizontal direction when the occupation probability is p. 
The function n v (p; L) is the probability, that the lattice percolates in the vertical direction, 
n hv (p; L) is the probability, that the lattice percolates in both directions, n s (p; L) is the 
probability, that the lattice percolates in the only direction (horizontal or vertical). It is 
clear, that on the square lattice TTh(p',L) = ir v (p;L) and 7Th(p',L) = \ix s {p;L) + ir hv (p;L). 
The example of functions ir s (p; L) for different L is shown on the Fig |T| a). 
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B. The (/-state Potts model as a correlated site-bond percolation. 



In the g-state Potts model the spin variable a takes values from the set {1, 2, ... , g}. The 
probability P(uj) of spin configuration to is defined by 

Ppotts(oo) = exp(-/3H(u))/Z(j3), Z(J3) = ^ Pp us(uj), H = -J ^ 5^ a ., (5) 

w <i,j> 

where H(u) - the Hamiltonian of the Potts model, Z{(3) is the partition function, and 
< i,j > is the sum over all neighbor sites. We put through this paper the coupling constant 
J = 1 and Boltzman factor fee = 1- 

Fortuin and Kastelleyn propose mapping of g-state Potts model onto the site-bond cor- 
related percolation JF/jl. It was used by Swendsen and Wang for their Monte-Carlo cluster 



algorithm jri| . The single-cluster algorithm was proposed by Wolff |19[ . Thermodynamical 
quantities can be expressed in terms of clusters in correlated site-bond percolation |2"0| , |2~T| ] . 

Lets the bond between neighbor sites with same values of spin variable Oi = Oj is closed 
bij = 1 with probability r = 1 — exp(— (3) and opened bij = with probability 1 — r, and 
the bond between neighbor sites with different values of spin variables is always opened. 

Define the chosen configuration of closed and opened bonds by v. Then joint distribution 
of the spin configuration u and the bond configuration v is (see p2]j21|) 



P Potts (uj,v) = Z- 1 (P) J] [(l-r^ + r^J, Z((3) = J2I2 P ( U} ^ ^ 

<i,j> u> v 

where Pp otts {^) = Ppottsiyi , v) is the probability of spin configuration u and Z[f3) 

is partition function for g-state Potts model. To each spin configuration u in accordance 
with (§) corresponds a bond configuration v, when bonds between neighbor sites with differ- 
ent values of spin variables are always opened, and the bonds between sites with equal values 
of spin variables are closed with probability r = 1 — exp(— /3) and opened with probability 
1 — r. The subset of sites, connected by closed bonds, called "physical" cluster. If we assign 
to spin variables in each "physical" cluster the random value from the set {1,2, ... ,q}, we 
obtain the new spin configuration. This procedure describes the Swendsen- Wang cluster 
algorithm, used in this paper. 



C. Crossing probabilities for the (/-state Potts model. 

For the percolation problem the crossing probabilities 7Th(p;L), ir v (p;L), 7r^„(p;L), 
7i s (p; L) are defined as indicator functions Ih(u), I v (u)) averaged over all site (or bond) 
configurations in accordance with ([!]), (fj), @, @ Crossing probabilities are probabilities 
for a system with linear size L to percolate according chosen rules at probability p of a site 
(or bond) to be occupied. 

For a g-state Potts model statistical weights of site-bond configuration (in terms of 
Fortuin-Kasteleyn mapping) are defined by expressions @, therefore it is natural to define 
the crossing probabilities through the indicator functions Ih(w,v), I v (u,v) averaged over all 
site and bond configurations. 
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FIG. 1. Crossing probabilities for percolation and Ising model. 
7T s {p', L) for percolation (a) and for the Ising model (b) 

n s (x) - normalized crossing probabilities for percolation (c) and for the Ising model (d) 
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■Kh(P;L) = y ^ Ppattsjv, v)Ih{v, v) 

U) v 



(7) 



where ?r^(/3; L) is a probability, that at least one percolating in the horizontal direction 
"physical" cluster exists, (3 = ■=, is the inverse temperature. By the same way others crossing 
probabilities are defined. 



III. SCALING OF CROSSING PROBABILITY FUNCTION ir s (P;L) FOR THE 
TWO-DIMENSIONAL Q-STATE POTTS MODEL. 

A. Features of approximation. 

We use the Swendsen-Wang algorithm to generate the different spin configurations. For 
each spin configuration we generate the bond configuration. Then we decompose the lattice 



into independent clusters of connected sites using Hoshen-Kopelman [23] algorithm. After 
that we analyze crossing properties Ih(uj,v), I v (lu,v) of this configuration. We average 
indicator functions Ih(u),v), I v (uj,v) over iV configurations. 

N N 

= j: E 7 >^)' *M = n S (1 " + ^ (1 " (8) 

i=i i=i 

We compute crossing probabilities for 40 — 50 values of (3 in the interval \j3 c — d(3(L),(3 c + 
d/3(L)], where the width of the interval was proportional d/3 ~ L~ y - see Fig. [I] b). The 
total number of configurations for every value of j3 are N = 2 x 10 4 — 15 x 10 4 . It should be 
noted, that (3 c (q = 2)= \og(V2 + 1) = 2f3 c (Ising) = 0.8813735870.... 



For percolation we generate site configurations by "grand canonical" method [gj] - every 
site is occupied with probability p and is empty with probability 1— p. For each concentration 
p we generate iV = 5 x 10 5 configurations. The number of analyzed configurations is the 
same for small and large lattice sizes. We use four-point shift-register random number 
generator with maximum length 9689 [ 25] . The crossing probabilities are not self-averaging 
by definition (||), and the variance of irh and n s not depend upon lattice size. 

Then we approximate crossing probabilities 7r s (/3; L) for different q by the Gauss function 

n s (L, P) ~ A(L) exp (-±(p - f3 c (L)) 2 B 2 (L)^ (9) 

Normalized crossing probabilities n s (x) looks like Gaussian - see. Fig. [I] c)-d). On this 
figures we placed normalized crossing probabilities n s (x) = ^2=^7f\ ' K s((f3 — /3 C (L))B(L); L) 
where x = ((3 — (3 C (L))B(L). For normalization we use parameters A(L), B(L), (3 C (L), which 
we obtained as a result of approximation of the numerical data by the expression (g). 

We dont approve, that the crossing probability is Gauss function, although they look 
similar Fig. |l] c), Fig. [I] d). We use this approximation like convenient method to compute 
the amplitude, the inverse variance and the maximum location of this function. 

As a result of approximations, for every investigated value of q and L we get the set of 
values A(L), B{L) (3 C {L). 
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B. Scaling of crossing probabilities. 
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FIG. 2. Amplitudes A(L;q) of crossing probabilities ir s ((3;L). 

We plot A as a function of L on Fig. |2|. For g < 4 we expect finite size corrections to limit 
values in the form aL~ x , as it was shown for percolation []26[ |. We found, that the functions 
A(L) for q = 2, 3, 4 could be well approximated by 



A(L) ~ A + al~", g < 4 



(10) 



Results of approximation plotted on Fig. |2] by lines. Second, third and fourth rows of the 
Table | represents result for amplitude A , and scaling terms a and x for q < 4. In the 
Table | below a row with values of q placed results of approximation for this q. The case 
q = 4 is intermediate of models with phase transitions of the first order q > 4 and the second 



order q < 4 [27 



We can expect, that for g > 4 the crossing probability tends to zero with L increasing, 
because the correlation length in the critical point £ c is finite [28]. Values of correlation 
length £ c for q > 4 in sixth row of Table | are taken from [28], We can also expect the 
different behavior of A(L) for cases L < £ c and L > £ c For q = 5,6, 8, 10 behavior of A(L) 
have another form - see Fig. 0. 



A{L) ~ Ao-alog(L), g>4, L < £ c 



(11) 



A(L) ~ A exp(-aL), g > 4, L > £ c (12) 

On this figure we use the logarithmic scale for L axis, and we see, that points for q = 5,6 
lay on straight lines. So we approximate data for q = 5, 6 by the logarithmic low (IT]). We 
put results of approximation into seventh and eighth rows of the Table |I[ Also we plote 
results of approximation on the Fig. |3] by lines. 
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FIG. 3. Amplitudes A(L;q) of crossing probabilities ir s (f3;L) in log(L) scale. 

But the data for q — 8,10 looks like straight lines, when we use the logarithmic scale 
for A(L) axis, while the L axis use the normal scale - see Fig. |j. Therefore we approximate 
A(L) for q — 8, 10 by the exponential low flT2|). Results of approximation for q = 8,10 by the 
exponential formula (|i~2D are placed in seventh and eights rows of Table | and are marked 
by stars *, and plotted on Fig. ^ Fig. f|by lines. 

We see, that for q = 5, 6 the correlation length is greater, then the linear size of our 
lattices L = 8, . . . , 256. But for the case q = 10 we examine region L > £ c . For q = 8 the 
value £ c = 23.8 and only two sizes L — 8, 16 are smaller, than £ c . 

We make additional calculations to check this approximation formulas (|TTD ([12]). We 
compute the crossing probability vr s (/5; L) for q — 5, 6, 8, 10 direct in critical points !3 c (q) = 
l°g(\/iq) + 1) H3- Results of this calculation presented on Fig. [| and Fig. ||. This figures 
corresponded to figures Fig. [| Fig. [| respectively 

The picture, which we can see in this figures, allowed us to make the following conclusion: 
for the case q > 4 the amplitude A(L) of crossing probability goes to zero by the logarithmic 
low ([Ll]) - see Fig. |3|, Fig. ||, when L < £ c goes by the exponential low (0) - see Fig. f|, 
Fig. [| when L > £ c . For the case q < 4 the amplitude A(L) goes to finite value by the 
power low (|T0"D. 

The expression for the inverse dispersion B in (§) follows from the scaling relation (for 
the case of the Gauss function (^|), parameter B is the inverse dispersion B = -p). The 
universality of the function 7r s (/?; L) implies, that the inverse dispersion is proportional the 
lattice size L in the power of y — - at least for cases of the second order phase transitions, 
when the critical index v of correlation length is defined - see formula (p^3|). Behavior of 
parameter B(L) as a function of linear size L is shown on the Fig. ^ in log-log scales. The 
points for different q lay on straight lines in accordance with (|I~3|) 

^- r B(L)^U (13) 
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FIG. 4. Amplitudes q) of crossing probabilities ir s ((3;L) in log(A) scale. 



In eleventh, twelfth, fourteenth, fifteenth rows Table | placed results for b and y of 
approximation of inverse width B(L) (0). In the tenth row we put analytical values of ^ 
for q < 4. We see, that for q < 4 the scaling parameter y is equal the inverse critical index 
v in agreement with scaling relations tt s (/3; L) = f(L»((3 — (3 C {L))). But for g = 4 we see, 
that y ^ I- For g > 4 we obtain values of thermal scaling exponents y, but it is unclear, 
how to relate it to another critical exponents. 

For the locations /3 C (L), maximum of the function ir s (/3;L), we expect deviations from 
limit values in the power form L~ z . For the percolation, the power z seems to be equal 
z = 1 + - 



1.75 



But the same deviation of maxima for magnetic susceptibility and 

specific heat of the Ising model is proportional to L» [ 29, 30 1 So we considered z as a free 
parameter and use for approximation formula ([14]). 

(3 C (L) ~ (3 c (oo) + cL- z (14) 

In seventeenth and twenty second rows of the Table | we put the precise values of critical 
points. For q = 2, 3, 4, 5, 6, 8, 10 we use expression f3 c (q) = log(.»/g+l). We use data from |26| 
for the critical point for percolation. 

The position of the maximum of the crossing probability TTp s (j3;L) on the lattice with 
linear size L goes to the critical point by power low (TA). We placed parameters of approx- 
imation (3 c (oo), c, z in eighteenth, nineteenth, twentieth, twenty-third, twenty-fourth and 
twenty-fifth rows in the Table |. 

There in good agreement between results of approximation and analytical values of criti- 
cal temperature. This values deviates only in fifth and fourth digits after the decimal point. 
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FIG. 5. Crossing probability at critical points ir s ((3 c ;L) in log(L) scale. 

IV. SUMMARY. 

• crossing probabilities for g-state Potts model are universal functions of scaling variable 
(f3 - p c )L$, looks like the Gauss function Aexp(-(/3 - (3 C (L)) 2 B 2 L^). 

• locations of maxima of crossing probabilities (3c(L; q) goes to phase transition points 
/9(oo; q) = \og(y/q + 1) when L goes to oo. 

• for the case q < 4 the thermal scaling index y — where z/ is the critical index of 
correlation length. For q = 4 the thermal index y ^ -. For g < 4 the limit of the 
crossing probability in the critical point is nonzero. 

• for q > 4 crossing probabilities in critical points go to zero by the logarithmic low, 
when L < £ c and by the exponential low, when L > £ c , where £ c is the correlation 
length in the critical point. 
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TABLE I. Results of approximation by (|lT])-(|lJ). 



1 


Q 


q = l 


q = 2 


q = 3 


q = A 


2 


A 


0.3538(3) 


0.3059(8) 


0.2610(12) 


0.17(5) 


3 


a 


-0.73(60) 


0.36(23) 


0.32(13) 


0.17(3) 


4 


X 


-2.1(4) 


1.48(33) 


1.20(20) 


0.3220(29) 


5 


q 


g = 5 


q = 6 


q = 8 


q = 10 


6 


e(/3 c ) - see. H 


2512.1 


158.9 


23.9 


10.6 


7 




0.277(3) 


0.301(15) 


0.198(27)* 


0.231(7)* 


8 


a 


0.0269(8) 


0.048(4) 


0.035(3)* 


0.076(1)* 


9 


q 


q = l 


q = 2 


q = 3 


g = 4 


10 


1. 


0.75 


1.0 


1.2 


1.5 


11 




0.7343(28) 


0.9988(19) 


1.1945(26) 


1.337(11) 


12 


6 


2.131(29) 


0.961(8) 


0.766(8) 


0.68(3) 


13 


q 


g = 5 


q = 6 


q = 8 


q = 10 


14 


2/ 


1.5284(85) 


1.654(33) 


1.865(13) 


1.886(19) 


15 


6 


0.485(16) 


0.450(45) 


0.306(13) 


0.360(18) 


16 


q 


q=l 


q = 2 


g = 3 


g = 4 


17 


/3 C exact 


p c = 0.592746 


0.881374** 


1.005053 


1.098612 


18 


/3 c (oo) 


Pc (oo) = 0.592731(28) 


0.881267(28) 


1.004990(19) 


1.098600(22) 


19 


c 


-0.333(12) 


-0.309(17) 


-0.527(24) 


-0.74(4) 


20 


z 


1.57(12) 


1.148(17) 


1.277(14) 


1.391(18) 


21 


q 


q = 5 


q = 6 


q = 8 


q= 10 


22 


/3 C exact 


1.174359 


1.238226 


1.342454 


1.426062 


23 


/3 c (oo) 


1.17440(6) 


1.23828(5) 


1.3427(5) 


1.42648(18) 


24 


c 


-1.114(24) 


-1.30(9) 


-1.5(8) 


-1.89(15) 


25 


z 


1.542(7) 


1.608(25) 


1.68(19) 


1.78(3) 



* - approximation by exponent - formula (|12j) 

* - p c ( q = 2) = log(V2 + 1) = 2f3(Ising) = 0.881373587... 
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